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Magnetic diffractive scattering of coherent metastable helium atomic beams from CoO(OOl) sur- 
faces reveal the presence of an anomalous enhancement in the intensity of the half-order antiferro- 
magnetic diffraction peak over the temperature range of 250-320K. Electronic cluster calculations 
based on Hartree-Fock states with configuration interactions, together with mean-field and Monte 
Carlo studies of this system suggest that this behavior is associated with low-lying surface spin 
excited states. 
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The magnetic properties of surfaces and thin films are 
currently a topic of much interest. This arises from the 
fact that these properties may be completely different 
from their bulk counterparts. For example, a pure sur- 
face magnetic ordering transition can occur at a tem- 
perature T* higher than the bulk transition tempera- 
ture Tg. A possible occurrence of such a phenomenon in 
gadolinium and terbium has been the subject of several 
papers 0,0 EI- More recently, our group reported on the 
critical behavior of antiferromagnetic (AFM) ordering on 
the NiO(OOf) surface^, as measured by metastable He 
beam (He* ) diffraction. These measurements have shown 
a surface Neel temperature higher than the bulk value. 
Further, the temperature dependence of the sublattice 
magnetization showed a cross-over behavior consistent 
with the extraordinary surface transition belonging to the 
universality class of the semi-infinite anisotropic Heisen- 
berg modelQ- 

In this Letter, we show that the AFM ordering on the 
CoO(OOI) surface, as measured by He* diffraction, ex- 
hibits an anomalous enhancement over the temperature 
range between 250 K and 320 K, straddling the bulk Neel 
temperature T^ of 290 K. At T^,, the enhancement is 
suppressed and a dip in the magnetic diffraction peak 
intensity is observed. 

CoO belongs to the family of AFM 3d transition-metal 
monoxides with a rocksalt structure in the paramag- 
netic phase. The spin-ordering in the bulk is character- 
ized by ferromagnetic (III) spin sheets that are stacked 
antiferromagnetically 6j . On the (001) surface, this ar- 
rangement leads to a (2x1) AFM spin structure. Our 
experimental He* diffraction measurements, shown in 
Fig. ^ display half-order magnetic peaks along the (10) 
and (01) directions, and no peak along the (11) 

direction, not shown, and support the existence of such 
an AFM structure on the CoO(001) surface. 

It has been established 0, @ that the lowest bulk 
crystal-field (CF) split orbital state, 4 T4 ( 4 Ti ff ), has an 
effective angular momentum I = 1 and is twelve-fold dc- 
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FIG. 1: Diffraction spectra measured along the (01) direction 
on CoO(OOl). 

generate. This degeneracy gives rise to a strong Jahn- 
Teller instability, which is stabilized via spin-orbit (SO) 
interactions, giving rise to a j = 1/2 ground-state (GS), 
and j — 3/2, 5/2 excited-state (ES) manifolds. In the GS, 
a hole occupies a Kf^-type orbital, which reduces the re- 
pulsion between a Co 2+ ion and its neighboring 2 ~ ion. 
This leads to an appreciable tetragonal contraction of 
~ 1.2%(at low temperatures)along one of the principal 
crystallographic axes[8j; about two orders of magnitude 
greater than the magnetostriction observed in the other 
monoxides 0. Recently we have carried out electronic 
cluster calculations employing Hartree-Fock states with 
configuration interactions for bulk CoO which confirm all 
the above features. A detailed description of the method- 
ology and results are given in R ef. Jl Oj. 

In earlier publications 0, Q, Il2l Il3l H3 | we demon- 
strated that He* diffraction is very sensitive to surface 
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AFM ordering. In particular, in Ref. |4| we have shown 
that the microscopic origin of He* magnetic scattering 
can be accounted for by an imaginary, spin-dependent 
component of the He*- surface interaction potential of 
the form AV s (r) oc i V[S(R);z] (S Z (R,T)), where 
(S Z (R, T)} is the thermally-averaged z-component of the 
spin at (R, z) and V [£(R); z] is the corresponding local 
surface density of states. Within the framework of the 
cikonal approximation, the diffractive scattering ampli- 
tude Ac, resulting from the above imaginary potential, 
together with a corrugated hard-wall real component of 
the surface potential, with a corrugation shape function 
C(R), is given by 

e -(W+a) r 

A G = / dn exp {i [G • R + q z ((R)] 

" Ju.c, 

-£(R) (S Z (R))} 

(1) 

where f2 is the area of the surface primitive mesh, W 
is the Debye- Waller factor and a is an attenuation fac- 
tor accounting for He* beam decay through the Penning 
ionization channel. £(R) includes P[5(R);z] plus other 
constants. In the case of the (2x1) spin-ordering, the 
diffraction amplitude of the i-order peak is related to 
the average sublattice magnetization by Q 

-A(iAO) = e-( w+a ^ h{x), (2) 

where x oc (S Z (R)) and Ix{x) is a modified Bessel func- 
tion of order one. 

The details of the design and operation of our exper- 
imental facility have been presented earlier in Ref. 
Beam intensities of about 3 x 10 5 He* atoms s _1 at the 
sample surface were achieved in the present measure- 
ments. All of the data were obtained from CoO(OOl) 
surfaces freshly cleaved in vacuum, with a background 
pressure less than 10 -10 Torr with He beam turned off. 
The CoO sample rods were attached to a copper sample 
holder with silver conducting epoxy. A silicon diode was 
attached at the base of the sample holder to measure 
and control the sample temperature. The temperature 
controller (Scientific Instruments Inc. Model 9600-1) is 
factory calibrated and is accurate to ±0.5 K. Scans were 
taken at 2 K increments between 250 K and 320 K and at 
larger increments at lower temperatures. Measurements 
of the diffraction intensities consistently yielded a mini- 
mum in the i-order magnetic diffraction peak intensity 
at T b N ~ 290K for all cleaved surfaces, irrespective of 
their distance from the silicon diode. We determined T^ 
to be 290±0.2 K by measuring the specific heat of several 
bulk samples as a function of temperature with a Pcrkin- 
Elmer model DSC7 differential scanning calorimeter. 

Fig-Elrepresents data collected from a single cleave and 
displays the enhancement of the peak intensity both be- 
low and above T^. Fig.[5]depicts data collected from sev- 
eral cleaves and shows more details of the enhancement 
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FIG. 2: Half-order peak intensity as a function of surface 
temperature. Data shown is from a single cleaved surface. 
The dashed line demarks the bulk Neel temperature. The 
solid line is a guide to the eye. 

below T b N . We notice that the peak intensity decreases 
slightly with increasing temperature from 50 K to about 
250 K, but then exhibits an anomalous increase with fur- 
ther increase in temperature, reaching a maximum at 
about 280 K, followed by a steep decrease, reaching a 
minimum at T^. Above T^ it increases once more to a 
maximum at 310 K and finally disappears at about 320 
K. This behavior has been checked repeatedly and was 
always reproducible even for different crystal ingots. 

To investigate these anomalous effects, first, we 
extended our electronic cluster calculations to the 
CoO(OOl) surface. Under the C± v symmetry of the sur- 
face, the state splits into a nondegenerate A2 and 
doubly-degenerate E orbitals, separated by 50 meV, com- 
parable to SO interaction energies, so that the latter can- 
not be treated as a perturbation on the CF state. A bal- 
anced calculation shows that all the ensuing states are 
Kramers doublets. The first ES doublets lie at 29 and 
94 meV above the GS. We also find that the expectation 
value of the spin in the GS and first ES to be is*) ~ 1/2 
and (s z ) ~ 3/2, with z normal to the surface. |10| 

Next, we constructed a model Hamiltonian that de- 
scribes the low-lying spin states and includes magne- 
toelastic effects. We started with the SO eigenstates 
of the Co +2 as a basis set and focused on the lowest- 
lying manifold 4 Fk/?- By using the operator-equivalent 
formalism (OEF)|l5j, a spin representation for the sur- 
face CF Hamiltonian is obtained as TLcf — cS z . Hcf 
splits the 4 F3/2 manifold into two doublets: rn s = ±3/2 
and rn s = ±1/2 with an energy gap of Ao = 2c = 29 
meV. We added the effect of the tetragonal lattice dis- 
tortion on the crystal field, by introducing a uniform sur- 
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FIG. 3: Intensity of the magnetic half-order diffraction peak 
as a function of surface temperature. Data shown is collected 
from several different cleaved surfaces, as indicated by solid 
and open circles and triangles. The dashed line demarks the 
bulk Neel temperature. The solid line is a guide to the eye. 



face tetragonal strain field rj, which, in the OEF, has 
the form bijS 2 , where b is the magnetoelastic constant, 
and modifies Hcf as Hcf — +br]) S 2 . The model 
Hamiltonian can then be written as 

ft = X^CF + J! ^CF + ^2 S Q • Sp 



<nnn> 



(3) 
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where a, (3 denote the two AFM sublattices and <j rep- 
resents a uniform bulk tetragonal stress field that cou- 
ples to the surface strain; it has the experimentally 
measured temperature-dependent form Ht| a = 
( t \ ' 6 

Go ( 1 — J . We evaluated the free energy of the sys- 
tem using the Bogoliubov-Peierls mean field variational 
method 18] with a trial noninteracting Hamiltonian of 
the form 



Hq — 22 ( H\ 
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where h a ^ are the effective mean fields along the pr- 
axis, taken as variational parameters. The temperature- 
dependence of the sublattice magnetization is then de- 
termined by solving the coupled self-consistent equations 
of the magnetization, M(T): 

_ 1 sinh(2/3J 2 M) + 3e~P A sinh(6/3 J 2 M) 
' ~ 2 cosh(2/3J 2 M ) + e^ A cosh(6/3J 2 M) ( ' 
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FIG. 4: Temperature-dependence of the half order peak in- 
tensity from a variational mean-field calculation. 



and of the gap, A(T) 
2ab 



A = A - 



b 2 cosh(2/3 J 2 M) + 9e- f3A cosh(6/3J 2 M) 
~ 2k cosh(2/3J 2 Af) + e-^ A cosh(6/3J 2 Af ) ' 

(6) 

These equations are obtained by minimizing the free en- 
ergy with respect to h ai hp and r\. The solution of 
these equations is carried out using experimental val- 
ues for a and k, as obtained from neutron scattering 
measurements [l9| of the phonon spectra, which yields 
k = 600A . 

The temperature-dependence of the half-order peak in- 
tensity I(i/2 ft){T) is determined by substituting the com- 
puted M(T) into Eq. 10, with W = 0. Fig. H shows 
1(1/2,0) (T) where both the anomalous enhancement and 
the depression at are manifested. It corresponds to 
J 2 /Ao = 0.25 and b = 8.8 An, comparable to reported 
experimental values for CoO |2fJ, |2l| . 

A Monte-Carlo simulation, based on the following 
Hamiltonian 

H = H s + H s f, + Hb 

H a = J a £ St-S^ + ^l^ + ba^lk) (S?J 2 

<nnn> i \ k / 

i k 
<rmn> 

(7) 

also confirms the anomalous effects. Here, the indices 
s and b refer to surface and bulk, respectively, and 
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FIG. 5: Computed surface magnetization (filled squares), 
bulk magnetization (open squares) and the excited state pop- 
ulation (filled circles) as a function of temperature. 



s z = 1/2 state. As the temperature increases, the bulk 
spin ordering decreases and the surface spins, less con- 
strained by the bulk, re-orient by populating the first 
s z = 3/2 excited state. The reorientation arises because 
of two contributions to the free energy: The energy of the 
3/2-state is lowered by an amount ~ S 2 J S through cou- 
pling to surface neighbors, with S = ±|, and the entropy 
is increased by populating new states. The depression in 
half-order peak intensity at Tjy is due to magnetoelastic 
coupling to the bulk. 

This work is supported by the U.S. Department of En- 
ergy under Grant No. DE-FG02-85ER45222. 



a = s or b. The sum over k runs over all the nearest- 
neighbor O -2 ions for a Co +2 ion at site i. We show typ- 
ical results in Fig.[S]for a 20 x 20 x 33 slab with periodic 
boundary conditions. Filled squares represent the mag- 
netic diffraction peak intensities for the surface, I S M {T) 
with the Debyc- Waller factor W = 0, while the open 
squares show the corresponding intensities for the bulk, 
I b M (T). Im(T) manifests both the anomalous enhance- 
ment and the dip at Tjy, similar to both experimental 
and mean field results. It also shows that > Tjy. 

The filled circles in Fig. 5 depict the population of the 
excited state n e , which mimics the behavior of /| f (T) 
in the anomalous region. n e is appreciably higher than 
would be expected from thermal activation across the gap 
Ao . We also note that the bulk spins do not populate the 
j = 3/2 manifold, since the energy gap is much higher 
than T%. 

Despite the fact that both mean field and Monte-Carlo 
results show a suppression in the intensity at T^-, a much 
more pronounced attenuation is observed in the experi- 
ment. This difference in the magnitude of attenuation 
can be attributed to the fact that the model we used 
is a simple one with the simplest choice for crystal field 
potential and magnetoelastic coupling. But it should be 
emphasized that even with the current model, the most 
important features of the experimental results, namely 
the enhancement and the supprcsion of the intensity at 
T^r are captured. 

In conclusion, the following picture emerges from the 
results presented above. At low enough temperatures, 
the ordered AFM bulk spins pin the surface spins to the 
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